NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



OPTIMAL DIGITAL CONTROL 
OF 

A 3ANK-T0-IUPII MISSILE 
by 

Carlos A. L. Velloso 
March 1984 

Thesis Advisor: Daniel J. Collins 

Approved for public release; distribution unlimited 



T21'570P 



U1 i^-;j_CLOO JLJ. _L^va 



SeCUUlTY CLASSIFICATION OF THIS PAGE 0»tm Enfrmd) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. PCPOPT NUMBCP 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG HUMBER 


4. TITLE rand SublllU) 

Optimal Digital Control of a Bank-to-Tum Missile 


5. TYPE OF REPORT A PERIOD COVERED 

Master Thesis 
March 1984 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHONCaJ 

Carlos Augusto Leal Velloso 


8. CONTRACT OR GRANT NUMBER^*; 


t. PENFONHINC OMCAnIZATION NAME AND AOONESS 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM element. PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


tl. CONTROLLING OFFICE NAME ANO AOORESS 

Naval Postgraduate School 
1 Monterey, California 93940 


12. REPORT DATE 

March 1984 


■ 3. number of paces ' 

154 , 


U. monitoring agency name a AOORESSrif dtitmrmnt trom Controiltng OHicm) 


IS. security class, (of thia report) 

Unclassified 


1S«. DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


14. OISTNiauTION STATEMENT (oi' lhl» Report) 

Approved for public release; distribution unlimited 


17. DISTRIBUTION STATEMENT (of tho oOotroet ontorod in Block 20, ii dlUoront trom Hoport) 


IS. supplementary notes 


IS* KEY WORDS (Continum on rororvo otdo ii noeoooorr tdontiiy by bioek numb9r) 

Guidance and Control 
Optinal Digital Control 
Terminal Guidance Law 
Bank-to-Tum Missile 


20. abstract (Conttnum on rmworom otdo ti nmcommmry ond idmnttty by btock numbmr) 

This work adresses the application of digital optimum control theory to 
a bank-to-turn missile. 

A optimal guidance law has been developed and tested in several scenarios 
using a 2-D model. Effects of sample rate, pitch angle, gravity and approxima- 
tions for small and large roll excursions are discussed. 



1 JAM*7I 1473 edition or 1 NOV •! IS OMOUETB 
S/N 0102- tF- 014- 6601 



Unclassified 



1 



SECUNITY CLASSIFICATION OF THIS PACE (Whan 0*>« inffd) 



Approved for public release; distribution unlimited. 



Optimal Digital Coctrol 
of 

a Bank-t o-Turn Missile 



by 



Carlos A.L. VelLoso 
Major, Braz ilian' Air Force 

B.S., Institute Tecnologico da Aaronautica, Brasil, 1976 



Submittad in parr ial -fulfillment of the 
reguirsmsnes for the dagrea of 



MASTER OF SCIENCE IN ENGINEERING SCIENCE 



f rom the 



NA7AL POSTGRADUATE SCHOOL 
March 1984 



r 

M/ 

v: 



&BSTHACT 

This work adrssses tha application of digi 
control theory to a bank-to-turn missila. 

A optimal guidance law has beea davelopad an 
ssvaral scenacios, using a 2-D modal. Effects 
rata, pitch angle, graviry and approximations fo 
large roll excursions are discussed. 



al optimum 

tested in 
of sample 
small and 



3 



TABLE OF CONTESTS 



I. INTRODUCTION 11 

II. MODEL OF THE SYSTEM 13 

A. INTRODUCTION 13 

3. ASSUMPTIONS 13 

C. THE CONTINUOUS MODEL 13 

D. THE DISCRETE MODEL 17 

1. Introductioa 17 

2. Calculctioa of the Matrices A(k), B(k) 

and E (k) 19 

III. THE OPTIMAL CONTROLLER 29 

A. DERIVATION OF THE OPTIMAL CONTROLLER 29 

3. EXTENTION OF THE MODEL FDR LARGE ROLL 

ANGLES 36 

1. Effects on the Miss Disrance of the 

Extention of che Model 4 2 

C. SIMULATION 44 

D. COMMENTS AND CONCLUSIONS 45 

1. Resales 45 

2. Connenes 47 

IV. ANALISYS OF GAINS, S AMPLE RATE AND PITCH ANGLE . . 82 

A. ANALISYS OF THE GAINS 82 

3. EFFECT OF THE SAMPLE RATE 94 

C. EFFECT OF THE INITIAL PUSH ANGLE 85 

D. EFFECT OF TIME TO INTERCEPT 88 

V. FINAL CONCLUSIONS AND COMMENTS 130 



4 



APPENDIX A: F03T3AN PROGRAM 132 

LIST OF REFERENCES 153 

INITIAL DISTRIBUTION LIST 154 



5 



. . t - , * . u j 

. * * - 'T?#4 ^^"Tr.iliy.'Ue 



TUIlEBtl 
I IfiMI 



LIST OF TABLES 



I. Results From Tesrs 81 

II. Rssulrs from Biased Control 96 

III. Results for Differants Sampis Paricds 105 

IV. Results Using Pitch Angle 120 

V. Effect of Time to Intercept 129 



o 



« ^ 




LIST OF FIGORES 



2. 1 Referencs Frames 27 

2.2 Refsrenca Frames 28 

3.1 Reprasenta t icn of the System 51 

3.2 Variation of commands wi-h Initial Roll 

Angle 52 

3.3 Variation of Roll Angle - Small Ang] es 53 

3.4 Variation of Roll Angle - Large Angles 54 

3.5 Ideal Initial Roll Angles 55 

3.6 Scenarios for Simulation 56 

3.7 Commanded Acceleration- Case 0 57 

3.8 Commanded Roll Rate- case D 53 

3.9 Miss Distance in Z Direction- Case 0 59 

3.10 Miss Distance in Y Direction- Case 0 60 

3.11 Roll Angle- Case 0 61 

3.12 Commanded Acceleration- Case 1 62 

3.13 Commanded Roll Rate- Case 1 63 

3.14 Miss Distance- Casel 54 

3.15 Roll Angle- Case 1 65 

3.16 Commanded Acceleration- Case 2 66 

3.17 Commanded Roll Rate- Case 2 67 

3.18 Miss Distance- Case 2 68 

3.19 Roll Angle- Case 2 69 

3.20 Commanded Acceleration- Case 3 70 

3.21 Commanded Roll Rate- Case 3 71 

3.22 Miss Distance- Casa3 72 

3.23 Roll Angle- Case 3 73 

3.24 Commanded Acceleration- Case 4 74 

3.25 Commanded Roll Rate- Case 4 75 



7 



3.26 Miss Distance- Case 4 76 

3.27 Roll Angle- Case 4 77 

3.28 Corrected Model 78 

3.29 Boundary of the Commanded Acceleration 79 

3.30 Corrections on the Roll Angle 80 

4.1 Gains affecting the commanded acceleration ... 90 

4.2 Gains Affecting the Commanded Roll Rate .... 91 

4.3 Commanded Acceleration-Case 5 92 

4.4 Commanded Roll Rate-Case 5 93 

4.5 Miss Distance-Case 5 94 

4.6 Roll Angle-Case 5 95 

4.7 Commanded Acceleration-Case 6 97 

4.8 Commanded Roll Rate-Case 6 98 

4.9 Miss Distance-Case 6 99 

4.10 Roll Angle-Case 6 100 

4.11 Commanded Acceleration-Case 7 101 

4.12 Commanded Roll Rate-Case 7 102 

4.13 Miss di sta nce-Case 7 103 

4.14 Roll Angle-Case 7 104 

4.15 Effect of Initial Pitch Angle 106 

4.16 Scenarios With Pitch Angle 107 

4.17 Commanded Acceleration-Case 8 108 

4.18 Commanded Roll Rate-Case 3 109 

4.19 Miss Distance-Case 8 110 

4.20 Roll Angle-Case 3 Ill 

4.21 Commanded Acceleration-Case 9 112 

4.22 Commanded Roll Rate-Case 9 113 

4.23 Miss Distance-Case 9 114 

4.24 Roll Angle-Case 9 115 

4.25 Commanded Acceleration-Case 10 116 

4.26 Commanded Roll Hate-Case 10 117 

4.27 Miss Distance-Case 10 118 

4.28 Roll Angle-Case 10 119 



3 



4.29 Commanded Acceleration-Case 11 121 

4.30 Commanded Hell Rate-Case 11 122 

4.31 Miss Distance-Case 11 123 

4.3 2 Roll Angle- Case 11 124 

4.33 Commanded Acceleration-Case 12 125 

4.34 Commanded Rcll Rate-Case 12 126 

4.35 Miss Distance-Case 12 127 

4.36 Roll Angle-Case 12 128 



9 



iCKNOSLEOGEMESr 



To tny wife. Rosy, for her constant, sncouragmsn t, supper 
and love. 

To my children. Ana, Antonio, Chica, Joao and Rita, fo 
yours smiles and love. 

Without you this worJc would be impossible. 



10 









tl’ 







...^ 

,r#f(iR 4 Jg 






.‘ Q*rsiratt -I iitrli « iiMf t^'ii'V t»* 




r»VOl bU 



•i% fee 

.*}Bn^ 



: .as >tviof«* .t«f .j» 36 U«'> rr .'t 



’ J 



,»fsl iff 1 * 1 .: Kaui>l 



.J 



. *i 3€J iii'aii •*-,■»% i *6r £i;v rnoxirl'fi 



lU 



■"V 





..y 



.... 









i. 





• * 



I. IN TR0D0C TI3M 



Because of threats from highly maneuverable high 

performance aircrafts and the need for increase standoff 
ranges, major improvements are needed in guidance and 
control capabilities of missiles. 

The high maneuverability of targets, has led to defense 
missiles capable of develop higher Lift accelerations and to 
more complex control laws, able to improve performance over 
well know laws as proportional navegation. 

In order to accomplish these new requirements with large 
standoff ranges, propulsion systems using airbreathing 
engines has been studied and developed in recent years. 

The advent of airbreathing engines leads natural to a 
consideration of bank-to-turn missiles in order to minimize 
the angle of attack of the inlets. 

The necessity of more complex control laws, leads in a 
general way, to the application of modern control and esti- 
mation theory, since more complete informations of the 
states of missile and target are necessary than those states 
informed by sensors commonly in use in missiles today. This 
leads to the use of a airborne computer. 

The present work adr esses the design and evaluation of a 
optimal digital control for application to terminal guidance 
in a bank-tc-turn missiles. 

One continuous two dimensional modal was adopted, in the 
following form: 



X (t)= A (t) X (t)+ B (t) u(t) ^ E g 



( 1 . 1 ) 



1 1 



1 

f 

whare the efface of gravity appears explicitly in the third 
term on the righ hand sida of axprassion 1.1. 

After the development of a equivalent discrete model, 
the optimal control problem has been solved, using a modified 
Ricatti equation due to the existence of the third term 
representing to the gravity effect. 

Next, several analysis has bean made in order to check 
the effects of small and large roll excursions, the affect 
of the sample rate on the system, and the effect of the 
initial pitch angle, in order to check the validity of such 
two dimensional model, when applied in some scenarios of 
interst. 
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II. MODEL OP THE STS^M 

A. INTfiODOCTIOH 

In "he present work the problem of terminal guidance for 
long range, bank-to-turn missiles with ramjet engines, using 
a digitalized system has been investigated. 

The model developed in reference 1 is used as the bass 
for this work. After the digitalization of that model, an 
optimal control law was developed. 

B. ASSOaPTIOHS 

Keeping the same assumptions as in ref. 1, on® has: 

The missil is limited to -2g*s and +15 g*s of commanded 
acceleration in the pitch plane, with zero lag. Also its yaw 
auto pilot has zero lag, yaw regulator maintains zero 
sideslip. 

Missile thrust exactly cancels drag. 

The angle of attack is assumed to be vary small, which 
leads to the commanded acceleration acting normal to the 
velocity vector. 

The missile will not have to roll through a large 
angl 9 . (Further considerations will give to this at the and 
of the derivation of the control law> . 

C. THE CONTIHHOOS MODEL 

Using the same reference frames as in ref. 1, one assumes: 
-Body frame with axis parallel to the longitudinal 
missile axis, positive y^^ axis out of the left wing ,and 
positive axis upward, (see fig. 2.1) 



-Fligth pa“h axis with Xp axis parallsl to th-a velocity 
vector , positivs Zp axis pointing upwards and yj. axis 
pointing to the left, (see fig. 2.2) 

In fig 2.1 and 2.2, the angles ^ and are the Eulerian 
roll and pitch angles. 

The state vector is given as 



i = [y. 



y. A^-y 



Atz 



A/] 



( 2 . 1 ) 



where y^ and Zp are the components of the relative target 
position, Yp and z, are the relative oarget velocity, Aty and 
Atz are the components of targen acceleration , which is expo- 
nentially decaying with a time constant <5 . 






( 2 . 2 ) 



where ^ is the initial roll angle (at t=0) . 
The control vector is given as; 




where Ac is the commanded acceleration and ?c is the 
commanded roll rate. 

The nonlinear plant equation is 



X =f ( X, u ) + G 



(2.4) 
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or 







0 


Aty i- Ac StN ^ 




0 


‘Aty/S 




0 


• 


•h 


o 


Atz - Ac cos ^ 




- ^ cot B 


-Atz /g 




0 


Pc 




0 



whsre g is gravity’s accslaration aad 0 is -h 
as seen in fig .2.2 

Linearizing and setting 

G = E g 
one has 





»• — 
• 

y p 




0 


0 




AtypAc fcotyl) 






0 




-Aty/S 




0 


0 


X - 


i. 


-f 


0 


0 




Atz > A'cCcos (f>) A(j> 




- ^^<5 ;4 
0 


o 




-Atz/s 




o 

0 




0 




0 


0 








J 



where 



( 2 . 5 ) 



e pitch angl 



( 2 . 6 ) 



( 2 . 7 ) 
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(2.3) 



- b 



0 0 0 -CCS 5 0 0 



1 

] 



in eqn. 2. 7, we have set 

cos^ = cos(^ =cos j6 cos2J)2^ - sin^ sLnAp 

sin ^ = sin(^ * ) =sin)^ cos Aj6*’Cosj6^ sin 

and expanded in which is considered small. 

Now assuming that A.^c, which is actually the desir 
control Ac, can be expressed in the form of: 






Ac = A CO 



>v] 



(2.9) 



with Aco = Ac at t=0, one has 



_x = AX'*’Bu + E g 



( 2 . 10 ) 



wher e 



A = 



/ 



0 0 



/ o 



0 0 o 



0 0 0 



O 0 o o 



0 o 
0 o 



00-^/3000 



/ 

0 

o 






Aco 



o 

/ Aco 



/- 



0 

t 

r; 

0 

O 



O 



o 



COi 



sm 



i 



( 2 . 11 ) 
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( 2 . 12 ) 



& 



0 0 

^>0 <f>Q 0 

0 0 

0 0 

- C05 0 

0 0 

0 1 



E 



0 0 3 0 



COS 6 




(2.13) 



g = g 



(2. 14) 



D. THE DISCRETE MODEL 
1 . Intro duction 

the inorcducricn of = Irgital ccmpurar no 
control the cont inuo us-time system, nne has to have some kind 
of intariacs in order to taka cara of the communication 
between the discrete and continuous- time systems. In this 
case it will be considered, A-to-D and D-to-A converters as 
samplers and zero-order holders as in reference 2. 

In such case, considering tae system; 

x(t) = A(t) x(t) + B(t) u(t) + S(t) g(t) (2.15) 
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one can vrite the state of the sysnem at t(!c + '; 



/ 






, tH*i 






where ^(t,to) is the transiction matrix of the syst 
ssntad by egn. 2.13. 

Furrher more , we will considered that the 
instants are equally spaced, or: 



t - t = r 

iC.I R 



K*> 



= kT + r 



so one can replace 
t = kT 

thus , 



x(kT+T)= )2^(kT+D x(kT) 



./cTfT 



+ M (kr^^T ) 3(^y|) ii(kr) 






/ k7*T 

^ (kr + T,-v| ) E(^ ) d^ g (kT) 



k-r 



or in a simplified notation: 



1 3 



(2. 16) 



sm repre- 
sa mpling 



( 2 . 17 ) 



(2.18) 



(2. 19) 



X (k + 1) = (k) x(k) + Bd (k) u (k) 



( 2 . 20 ) 



+ E^(k) g(k) 



wh a r e , 



A4 (k) = p (kr + T,kD 




KT^^r 



Bu (k) = 0 (kT+T, /n ) B(<n ) d 



"1 

IfT 



T| ) a-i^ 



^KTtT 

-d (kT+r,^ ) B(^ ) 

kT 



2. Calculation ths Mat^cas A (k) , 3(k) 



It is straightforward to show, using ti 
of the matrix A (t) , that the transicticn matri; 
2-19 is : 



<^(AThT) = 



/ 


r 


Ad 


0 


0 


0 




0 


{ 


'^^,3 


0 


0 


0 


A<Jz.t 


0 


0 




0 


0 


0 


0 


0 


0 


0 


f 


T 






0 


0 


0 


0 


/ 






0 


0 


0 


0 


c? 




0 


0 


0 


0 


0 


0 


0 


/ 



(2.21) 



( 2 . 22 ) 



(2.23) 



nd E (k) 

a sparseness 
of aquation 



(2. 24) 
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wher e : 

Ad. , =Ad =5 I - 5 2 (1-e'^ ) 






Ac) =Ad = (1-e'^C) 

'/.J '' 5.6 



For the calculation of the others terms one may mak 
use of the property of the transiction matrix that: 



dt. 



= A(t_2 ) Ad(t^ ,t, ) 



so. 



d Aj (kT+7, icT ) 

= A(kT+T) Ad (kT+T,kr» 

d^kTi-T) 



For A (2 ,7) : 



dCkT^T) 



= j (kT+T) = Aco 



, . A. 

Ti 






kr^r 



^ X CO 

6(kTu7) 

■ KT*T 



1- tZiI- 

~~T~~ 



:03 






Aj^^= Aco cos^ j ^1- i(kT^T) 



kr 



Ad = Aco CDS 
r 



A ["-[-yrr-] 



For A (1,7) : 

d Ad,., (kT+T, kl) _ 



d rkT+r) 



Ad^^(kT+T, kT) 



KT*T 



Ad (kT+T,kT) = I Ad d(kT + T) 

/ zj 

kT 
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T3 


j Aco COS^ 








Doing th<= same 


process for 


Ad (5, 


7) and 


Ad (U ,7) one 


has: 














Aco sin;Z< j 


12- 




i T3j 




A = Aco sin^ j 


I 


-1 


■ 2 ki- 1 " 

L . 


j T2 


1 


For nha ds- 


rivation of the 


rna' 


tr 


ix Bj (k) one needs 


according to egn. 2 


-22 












(kT+T) B( -j ) = 


Ad 


(kT+r,')| ) B(ij ) 






wher e. 
















* / 


At, Aft, 0 i 


0 


0 


*i,7 








/ 


A/n 0 

'Z.i 


0 


0 


Aff\ 

*2.r 






0 


0 


Ann 0 

h.i 


0 


0 


0 




A^(kUT,^ ) = 


0 


0 


0 / . 


A-tj 

U.s 


A 01 Arf, 




0 


0 


0 0 


/ 








0 


0 


0 0 


0 


A 


0 

<> 






0 

. 


0 


0 0 


0 




0 1 




where A<^ represents 




(kT+T ,-yj), and 








0 

0 

0 

0 

0 

0 

0 
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thus 










5in 




0 


0 






-cos 


A»n 

<Sr 


0 


0 


0 


0 



whsre 



A/»j 



ul 



= kT>T- ^ 5 



A 



2.T 



.m 

M.r 



I kT+r- i-ItZl— - 1 

L 2T\ < 2 t; J 



Ago cos 



= [t<w»ii - 



(kr^D T2 - 



T 4fi + 


2k ^• 1 


T2 /ri 


1 


1 


2 Ti 


J 1 


1 J 



Ago gos 






^'Tl = 
< 4,7 



• 


' 2k + 1 


T(kr + T) - 


Z Ti 



(kH-I) T2 - 



- T -T) <• 



• "■ 
2k + < 


- 


2t; 


T2 /Y| 



AGO 






Am = 
' 5,7 



(kTrT) ^ -,2 

kT+r- ~ ^ + -j — 

2 Ti ' ^ T.- 



Ago sin 



and for B^(k): 
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B, 



wher e 



kTtr 



Ba 






Bd 




Jin 




0 


0 


Ba, 

4,1 




- Coi ^ 


Bd 

s,t 


0 


0 


0 


T 


r- 'Tj ) doTj 


sin = 



(2. 25) 



KT 



T‘ 

/ kT^T 

H.rj 



which after some algebric work has been found 



3 (k) = 



■l-i 



T" 


(kTi. t)^ 


(kr+r/ kT 


/kt)^ I 


L 


JTi 


r.- 


6 TI 



Acc cos ^ 



and 



.kT* T 






kT 



which can be found to be: 
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thus 









■ 


1 




X4 


7 ’ 




2T\ 

m mJ 


’ .il 


1 

... 


2ki. \ 


T4 


I 







Aco cos 



A CO sin ^ 



(kT + j) fxT+t)*. (kT)* 

....... ....^... kT - ---- 



A CO sin^ 



In the same way: 

0 

0 

0 



^(kTi.r,ffi) £(^) - 



-T C06 9 

- cos 9 



0 

O 



thus E (k) is equal to 

E=^0 00 0 - -j- cos 9 -r cos 6) 0 0 

where we ha^e considered ^ as a constant angle. (Further 

comments cn this after development of the control law). 

Notice that througout this work ,the commanded 

acceleration has been considered an unknown and the assump- 

tion has been made that it will be in the form of eqn.2.9. 
One will need in further developments to consider the 
control Ac as a known Ac (k) , which will be a constant between 
kT and kT+T. With such assumptions, the discrete represen- 
tation of the system is easily found to be : 
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X (k + 1 ) = A (k) X (k) + E (k) u (k) E g 



(2. 26) 



whare 







/ 


T 


Ad 


0 


0 


0 








0 


/ 


>4cj 


0 


0 


0 








0 




-Vg 

c 


0 


0 


0 


0 


= 




0 


(P 


0 


i 


T 










0 


C> 


0 


0 


1 




T 






0 


0 


0 


0 


0 


-/6 


’ 0 






0 


0 


0 


0 


0 


0 


i 


with 




















3 T- 


' 32 


(1- 


e-'^/G 


) = 










3(1- s 




= Ad 












-l! 

’i” 


Ac 


cos ^ 


4 












T Ac CDS ^ 














.li 

z 


Ac 


sin^ 


4 











(2.27) 



Ad^j= T AC sin 
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and 



Bj(i<) : 



and 



£d Ck) = 



5^ 

' o 



T 

0 

- JL. cos ^ 
Z 

- T cos 

0 

0 



0 

0 

0 

_ Xf cos © 
2 

- T cos 0 
0 
0 



t3 y 

-1- Ac cos ® 

6 



— Ac cos _ 

2 AO 



-]“3 y 

— Ac 5in^^ 



— Ac sm 



0 



T 
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Figure 2.1 Referenrs Frames. 
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Figure 2.2 Refereare Frames. 
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III. OPTIMAL CONTROLLER 

A. DERIVATION OF THE OPTIMAL CONTROLLER 

In order to have a suitable guidance law zo 
the control commands, we will niaimizing the 
performance index: 



H.l 



J= -j- x''(n)H (n)x (n) + ^ u'(k) 3 ( )c) u (Ic) 







where x (N) is the final state at t = Ti 

As we want to minimize the final miss dis 
weighting matrix W(N) is taken as 



/ 

0 

0 

0 

0 

o 

0 



o 

o 

0 

o 

0 

0 

0 



0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0/000 

o 0 O 0 o 

o o O 0 o 

0 O 0 O 0 



and Q(k) is a two by two positive definite 
weighting matrix to be chosen. 

In the derivation of the soluzion, reference 3 
followed keeping in mind zhat the szaze equation 
form : 



i mplement 
following 



( 3 . 1 ) 



z ance, the 



( 3 . 2 ) 



symmetric 

has been 
has zh^ 
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(3.3) 



x(k + 1)= A (M X (k) + B(k) u(k) + B g 



or 



X (k + 1) = f (X (k) ,u (k) ,g) 



(3.4) 



whsre rhe third term, which rapresaars rha affect of th< 
gravity has baan considered constant. 

Considering that the parfcrmanca index is in the form: 



J = J2^[x (N) J +\ L (k) I x(k) ,u (k) ,g(k)J 



(3.5) 



* = o 



we need to find a saguenca of u(k) that minimizes J. 

Adjoin tha system equation to J with a multiplier (n) 



M- I 



r= ^ l^x (J 



(N) 



+) \ L (k) X (k) ,u (k) ,g 



f r 


] 1 


1 fj^ lx (k) ,u (k) ,g 


- x(kti) 1 



+ A (k+1) I f^ 

and defining a scalar seguanca H (k) 
H (k)= L(k) X (k) ,u (k) ,g j + 

+ A(k+1) f^ X (k) ,u ( k) ,g J 

k=0, 1 ,2, n-l 



(3.6) 



(3.7) 
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one has; 



J= /^l^x(N) j - A\n) x(N) 



/v-» 




+ ) H(5c) - A(1C) x(k) + H(0) 



J 



Ks I 



Considering differencial changes in J: 



(3.8) 



dJ = 



[ ™] 



dx(N) + 



(3.9) 



N-l 



k= I 






^o(k) 



+ du(0) 

^ u( 0 ) 



cJHfo) 



+ du(k) ) * dx ( k) -f 



^ x(o) 



choosing the multiplier A(k) so that 



A'(k) 



=0 

x(i«r^ 



(3. 10) 



thus 



+ A^(kfl) -ti. = 

c) x(l<) <ix(lc) 
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and 



e) u(lc) 



= 0 



( 3 . 11 ) 



with boundary condition 



X(N) 



= M.. 

c) x.^w) 



(3.12) 



we obtain the minimization of the parforaanca index. 
In the oresent case we have: 






J= ~ x^(N) W(ll) X (N) > ^ u^ (k) 2(K) a(k) 

+ A^(k+1 ) A (k) X (k) +B (k) u (k) +E 3-x(k + 1) j 



(3. 13) 



and H (k) : 



H (k) = Y u"^ (k) 2 (k) u(k) * 



(3. 14) 



+ X (k^-1) 



A (k) x(k) + 3(k) a(k) + E g 



then in order to minimize J: 



(k) Q (k) + X'(k<-1) 3(k) =0 



(3. 15) 
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T 





or considering that Q (k) =Q ( k) 

Q (k) u (k) = - (k) X (k* 1) 

and 

)^(k) = (3.16) 

X Ck) 

>!(k) = )J^(k + 1) A(k) 

and 

^(N) = x'(N) W (N) (3, 17) 



Notice rhat we are not weighting the snares in rhe 
performance index, except the las- state. A more general form 
could be ob*ained, with all stares being weigrhing , if we 
change the eqn 3.16 to: 



>[(k) = X'(k+1) A(k) + x^(k) W1(k) (3.13) 

where W1 (k) is the weighting matrix of the states. 

With equations 3.15,3.16 and 3.17 one is able to find 
the sequence of u(k) that will give the minima controls. 

Such set of aquatics can be solved by the sweep method 
as in ref. 2 . 

We will look for a solution of the form: 



u(k)=-F(k)x(k) -FG (k)g(k) 



(3.19) 
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what means that the commanded acceleration and roll rate, 
will have a correction due to the sffecr of gravity. 

Placing : 



X(k) = S(k) i(k) + SG(k) g(k) 



(3.20) 



from eqn. 3.15 

Q (k) u (k) =-B^ (k) js (k+1 ) x(k+ 1) +3 3 (k + 1) g (k + 1> j 
from eqn. 3.3 



Q(k) u (k) = -B (k) S(k+1) A (k) x(k)i- 



[ 



(3. 21) 



+B(k) u(k) + S gj - B^(k) SG(k + 1) g(k + 1) 



so 



Q (k) + B (k) S(k+ 1) 3 (k) 



u(k) = 



-B^ (k) S(k+1) A(k) x(k) - B^ (k) S(k + 1) E g (k) 



-b"^ (k) SG (k+ 1) g (k + 1) 



considering that g is a constant 



u(k) =-[o(k) + B^(k) S(k+1) B (k)j 



(3.22) 
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+ 



(k) s (k+ 1 ) A (k) 
J (k) S (k+ 1) E 



X(K) V 

+ B^ (k) Sg(k+1) 



g 



] 



so ; 

u (k) =-F (k) x(k) -FG (k) g 

whare 



F 



[Q(k) + 3 ' (k) S (k+1) B(k) 




(3.23) 



• B^(k) s (k+1) A (k) 



FG 



= [ 



Q(k) 



(k) 



S(k+1) B(k) 



'J 



(3.24) 



[3^(k) S(k + 1) S + B^Sg(k+1)j 
from eqn.3. 16 and 3.19 

X(k)=A^(k) X(k+1)=A^(k) S (k+ 1) x(k + 1) +SG (k + 1) g (k+ 1)j = 

=A^ (k) S(k + 1) A(k) X (k) + A^ (k) S(k+1) B(k) u(k) 

+ a"^ (k) S(k+1) Eg + A^^(k) SG(k+1) g(k+1) 

from eqn . 3.19 

X(k) = A^(k) S(k + 1) A(k) x(k) + A^ (k) S(k + 1) 3 (k) . 
*[-F(k) x(k) - FG(k) gj + A^(k) S(k+1) Eg + 

+ A SG (k) g (k+1) 
so, as g is a constant: 

S (k) X (k) +SG (k) g= j^A^ (k) S (k+1) A (k) -A (k) S (k+1 )B (k) F(k)'j 
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•thus 



(k) S(k+1) E - A^(k+1) 3(k) FG(k) 

+ A (k) SG (k + 1) ] g 



S(k)=A^ (k)S(k+1) A (k)-A (k) S (k+ 1) 3(k) F (k) 



(3 



and 



SG (k) =A^ (k)S (k + 1) E-A^ (k) S (k+ 1 ) 3 (k) F G ( k) + (3 

+ A^ (k) SG(k+ 1) 

These equations, 3.25, 3.26, 3.23 and 3.24 can be 

backwards with the final condition: 

S (N) =W (N) 

SG (N) =0 

Notice that this satisfies our previous boundary 
tion in eqn. 3.17 where: 

(N) =W (N) X (N) 

S (N) X (N) =W (N) X (N) 



so 



S (N) =W (N) 

B. EXTEHTION OF THE MODEL FOR LARGE ROLL ANGLES 

Op to this point, one has to take into accoun 
throughout the development of this work, the angle ha 
considered small, it is necessary to relax this restri 



. 25) 



. 26) 



solved 



con di- 



t that 
s been 
ction. 
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In order to do that, the system has been broken in two 
blocks as in figure 3.1. 

The first block is a representation of the algorithm 
which will calculate the optimal commands. The algorithm has 
contain with itself an exact model of the system or missile. 
The model of the missile is initialized from the information 
on the initial states, the initial input command Aco and 
initial roll angle (at t=0 ) ; and computes the optimal gains 
and further the optimal commands which will be feed to the 
miss ile. 

The method adopted in computing the optimal commands is 
more easily understood if one considers figure 3.2. 

In figure 3.2, the lines numbered as 0 in the graph for 
Ac and for Pc are the optimal commands for a given initial 
roll angle ♦ Lines number 1 are the commands for a second 
initial roll angle(;^, ) larger than^, and so on. Thus in 
figure 3.2 one has a family of optimal commands for any 
initial roll angle. 

Notice that the upper line of the graph of Ac represents 
the accelerations of a missile whioh had at t=0 a correct 
initial roll angle in order to hit the target with no 
comands in roll rate. 

In the present method the computer performs the calcula- 
tion of the commands only for the first step of time and 
then feeds these commands to the missile. The missile is 
then driven to the next state (x(k-*-1|) and feed-backs to the 
computer the information on the roil angle at that step. The 
roll angle feed-backs from the missile is considered by the 
algoritm as the initial roll angle at t=0 and the next 
commands are calculated. Not ice that at this second step the 
algorithm will feed to the system the second command (at 
t=t, ) . This process is them repeated until t is equal to 
the intercepted time. 



37 



It is important to realize that with -his merhcd of 
calculation, since the algorithm was developed wirh the 
assumption of small roll excursions some error is expected 
due to the fact that in computing the gains by solving a 
Bicatti equation backwards, as has been done, it is necessary 
to update the system from t=0 to T.= Ii at each seep, and in 
this process the roll angle is not small. Notice however that 
we are applying the commands only in one step, and if one 
expects that the roll rate will decrease to zero, as we are 
increasing in time, the variation of the roll angle will 
tend to decrease, sc , we can expecned than the error will 
decrease as the time increases. 

Anorher important point tc be studied is how rhe missile 
itself (second block in fig. 3. 2 ) has zo be implemented in 
order to be valid for small and large roll excursions. 

It is considered chat one has the perfect knowledge of 
the commands, thus the missile is modeled as a stats vari- 
able system as in sgn 2.26 with the initial roll angle being 
update at each step. In nhis way the system will rake ohs 
initial roll angle as nhe summation of all previous ininial 
roll angles, (see fig 3.U) 

From the original variables one has for the stare , 
the following: 



A^(k<-1) =<d/(k)+ ?c(k) T 



(3. 27) 



which is shew in figg- 3.3, where the initial roll angle is 
kept constant and is update each step. 

Notice however that considering large roll excursions 
(see fig. 3.4) and keeping in mind that the angle has 

bean defined as 
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the expression 3.27 is not valid, since in acdeling the 



This leads to a change in the expression for the varia- 
tion of as in fig. 3.4. In fig. 3.4, the angle is 

update at each step, so one has: 



By consideration of the equation 2.27 it can be seen 
that by setting the element A (7, 7) to zero, one can obtaim 
equation 3.28. 

Thus the missile model for large roll excursions vill be 
represented by the following equation. 

x(k + 1)= A (k) X (k) + B(k) u (k) + S g 



system it was assumed that would be samll 




Pc(k) T 



(3. 28) 



where 



I T A 



0 0 0 A,, 



(3.29) 




O O A 



0 O 



O 



AO() = 0 0 O / T 



0 O 0 O / A^^ 

0 0 0 O O O 



0 0 0 0 0 0 
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whsre A(7,7)=0 ,for large roll excursions 
&(7,7)=1 ,for small roll excursions 





X' 

6 


Ac cos ^ 


r 5.0 


Tf 

1 


Ac cos^ 


0 




0 


- Jlf C 03 ^ 


Tf 


Ac 


* T cos <5 


SI 

Z. 


Ac 5'n^ 


0 




0 


0 




T 



(3. 30) 



0 

0 

0 

.If cos e 

z 

- r cos 0 
0 
0 



(3. 31) 



We will redefine rhe states x, , as rhe relative posi- 
tion, Xt as rhe relative velocity in the Y direction and x, 
as the rarget acceleration in Y direcrion. The stares x^ , 
x^ and Xj. has rhe same meaning, but in the Z direction and 
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is . In this representation of the system, x, (k) i 
equal to the component of the miss distance along the Y 
direction and: 



X, (k+1)= Y^(k+1)= y^(k)+ y^ (k) T + (3.32) 

+ 1^5 r- 52 (1- e )J Aty(k) Ac(k) cos^^ ^/^(k) 

T*- / T* V 

+ -j- sin^ Ac(k)■^ Ac(ki cosp^ Pc(k) 

The first three terms in the RH3 of eqn.3.32, are easily 

seen as the contribuition to the miss distance of respec- 
tively rhe previous miss distance, relative velocity and 
target acceleration. The following two terms represents the 
cont ribuiticn of the commanded acceleration and the last 
term represents the effect of coupled Ac and Pc, and tends 
to be small due to the cube cf the sample period. 

For the component of miss distance in the Z direction, 
one has: 

x^ (k-H)= Zp(k+1)= z.(k)<- Zp (k) T-^ 

+ [s T- 3 2 (1 - e’^ ) Atz (k) -j- Ac(k) A^{'<) - 

- -j- cos^ Ac(k)+ Ac(ki sin^ Pc (k) - -j- ccs 0 g 

where its terms have the same physical meaning as in the 
expression for x (k^-1), with the effect of the gravity added 
to the expression. 

Since one can notice that in tae representation of the 
miss distance in T direction appear two terms as a function 
cf cos , and in the rep resentaction of the miss distance 
in Z direction appear two terms as function of sin^ , it is 



intaresring to verify that the foarth term in the RHS of 
beta expressions acts like a correction for the fifey csrm, 
Eefering to fig. 3.4b, one can see that at any step of rime, 
the commanded acceleration is actually. 



The same idea can be applied to the expression for 

X (k+ 1) . 

The terms x^ and Xy , represent the relative velocity, 
and are: 



Where the two first terms in the RHS represents the effecc 
of rhe velocity and acceleration at a previous step, and the 
other rhree terms has the same meaning as previously stated. 

The terms Xj and x^ are the target accelerations, in 
this model being exponentially decaying. 

^ • E ffec ts on the Miss Dist anc a of the Sxten ti on of the 
Model 

In previous subsection, a extention of the model for 
large roll excurssions has been performed. Notice that there 
are two models of the system being used. The first one, used 




and considering small angles: 
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in the algorithm is valid only for samll roll excursions, 
and a second model, valid for small and large roll excur- 
sions used as a representation of the missile. 

The algorithm with the ficsr model, as explained 
before, is initialized at each step with the acrual roll 
angle of the missil, and performs rhe calculation of the 
comands. 

In order to check the effect of the extention of the 
model on the miss distance, one can define a ideal initial 
roll angle ( ), as the roll angle at t=0 in order to 

have the comanded acceleration vector pointing to the 
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roll angle after some time, and due to the vertical compo- 
nent of the accelsr ation, when this occurs the missile would 
be in a position above the ideal trajectory, which means 
that it would follow a course parallel to the ideal trajec- 
tory to intercept. 
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One can see than such problem will lead to a large 
miss distance, in the case that the narget acceleranicn it 
is not small. Thus, some correction is necessary in order to 
improve the missile performance. 

Figure 3.5b shows the missil at some point of its 
trajectory where it has reach its maximum roll angle, nhus 
it is at a parallel course with its ideal trajectory to 
intercept. At this point, since all the states are known, 

in is possible to recompute a new ^ ,^ea\ * this 

point the computer is feed with winh the states an this 
point it will compute the commands in order to drive the 
missil to the new ^,aea\ » which will innroduce the desired 
correction. Notice that such correcnion can be made during 
all the flight, from t=0 , to t=Ti. 

In the present work this method has been accom- 
plished by feeding-back to the computer the roll angle and 
all the states of the missile. With this information the 
computer is able to perform the calculation of the corrected 
commands at each step of time. howeveras the states are 
being updated, it is also necessary to update the time to 
intercept, which has been done using the time to go to 
intercept, or: 

Ti (k) = Ti- k T 

where Ti is the nominal time to intercept. 

C. SlflOLATIOH 

In order to keep the same assumptions as reference 1, 
the matrix Q (weighting the control) has been put as 
suggested in such reference or: 
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o] 



(3.33) 



'b, 

0 

with 

b, =5.78 103 

=5. 0 meters 

Five different cases were run: 

Case 0, tested with one simple model valid only for 
small roll angles begins with missil and target on parallel 
courses to the inertial x axis, the target 100 meters above 
the missile and with an evasive manoever exponentially 
decaying with time constant of 20 seconds. The initial 
acceleration of target was -.5 g's in y direction. (see fig. 
3.6) 

Case 1, with the same scenarios as case 0, but was run 
using the algorithm for large roll excursions. 

Case 2, is the same scenario as in case 1, except that 
the initial acceleration of target was -1.0 g's in the y 
direction. (Same as case 1 in ref. 1| . 

Case 3, the same as case2,with target acceleration of -U 
g's. ( Same as case 2 in the reference 1) 

Case 4, same as case 3 but with target at initial posi- 
tion 600 meters bellow the missile. (S ame as case 3 in refer- 
ence 1.) 



D. COMMENTS AND CONCLOSIONS 
1 . Results 



In 


case 0, 


the missile 


begins its trajectory 


commanding 


26.5 m/sec2 


and the time 


history of the control 


Ac follows 


exactly a straight line 


in the form suggested in 
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ref. 1, as seen in fig.3. 7The control Pc, baginir/g at . 
rad per second, is decaying and ceachs zero at k=80 {see 
fig. 3. 8). Figares 3.9 and 3.10 show the miss distance, where 
one can see that the missile is crossing the target with a 
CG-to-CG distance of 1.5 meters. 





Figured. 


11, shows the 


time history for 


the roll 


angle. 


which as 


expect, reachs 


a constant value. 


with the 


missile 


cro ssing 


target at t=Ti, 


with a bank of .44 


radia ns. 



In case 1, which was run with the model for large 
roll excursions, the missile kept the same Aco, but there is 
a very small increase in further commanded acceleration in 
order to correct the effect of the roll angle on the 
vertical component of the control Ac, as seen in fig. 3. 12, 
and table I. 

Figure 3.13 shows that tie roll rate decreases 
almost as before, and the final roll angle is .42 rds, as 
shown in fig.3. 15 and table I. Refecing to figure 3.14 rhere 
is a change in the final miss distance, which is better than 
case 0, due to a improvement in its Z component (see table 
I) . This results in a final CG-to-CG miss distance of .65 



mere rs. 

In case 2, the missile has the same Ac at 0=0, but 
with the correction for roll angle being increased due to 
the increase of target accelerarion (see fig.3. 16), the 
commanded Ac reachs a larger peak value. 

The initial value of the commanded roll rats is .69 
radians, which is larger rhan case 1, due to the increase in 
the target acceleration. In figure 3.13, chere is no 
noticeable change in the shape of me curve for Z direction, 
and in nhe Y direction the final disrance is about the same 
as in case 1. The CG-to-CG distance at r=Ti is .73 merers. 
The larger roll rats leads to larger roll angles as seen in 



4:4 



ig. 3.19, where the final roll angle is .72 radians. 
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The effect of target acceleration can be easily seer, 
in case 3, where one can see that with the same kco, the 
accelerations are largslly increased from this point, and 

the missil begins its trajectory with very high roll 

rate (see figures 3.20 and3. 21) . There is a change in the Z 
component of the miss distance, that decreases its final 
value to ,05 meters, but now the miss distance in Y direc- 
tion is made worse as shown in figure 3.22, which leads to a 
final CG-to-C3 distance of 1.5 meters. In fig. 3. 23 one can 
sea that the missile crosses the target with a bank angle of 
1.26 radians. 

In case 4, due to the position of target 600 meters 
under the missile, the initial commanded acceleration is 

negative and reachs the limit of -2 g's. The inirial roll 

rate begins at a smaller value than in case 3 but increases 
during the inirial part of the flight reaching its peak 
value at 1.75 seconds when again as in the previous cases 
begins to decay. As the missile banks tc roll angles larger 
than 90 degrees, the acceleration goes to positive values, 
as seen in figures 3.24 and 3.25 Figure 3.26, shows the 
worse case among these in respect to the miss disrance, 
mainly in the Z direction, and in the final cg-to-cg 
distance, which is equal to 4.43 necers. Also the final 
roll angle of 3.0 radians is the Larges* among all these 
cases, as seen in fig. 3.27. 

2 . Co m me n t s 

Defining the projected zero efforr miss distance 
(ZSM) as the miss disrance the missile cross the target with 
no commands. It can be calculate! at t=0 as the initial 
distance between target and missile plus the miss distance 
due to the gravity or: 
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ZEM 



z (0) + 



g t2 



(3.3U) 



I 



which is equal in all the three first cases. 

In the first three cases, the initial missile’s 
commanded acceleration is the same as seen in table 3.1. 
Considering that the control Ac necessary to correct the 
initial miss distance can be calculated as: 



Ac 






ze M 





26.7 



(3. 35) 



which is close to the initial control Ac. 

This suggests that the initial Ac, would be that cne 
necessary tc correct the initial ZEM in Z direction, which 
agrees with reference 1. Notice however that in all cases 
the initial Ac is less than the calculated value of Aco, 
which agrees with the previous statement that some error was 
expected in the initial part of the computations. 

Considering the ideal initial roll angle as defined 
before, one has; 






ideai 



tan'i 



Z£ 



(3.36) 



From table I, and figures showing roll angles, it 
can be verified that the missile is banking to reach angles 
larger than r in order to correct its trajectory to 

hit the target. 

Therefore, in all cases, the mrssile begins its 
trajectory with an Aco (discussed before) , and a roil rate 
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which is proportional to the targst acceleration in 
direction. As the missile rolls at decreasing roll rates, 
the ccmmanded acceleration is changed in order for co 
compensate the effect of the roll angle on its Z direction 
component. At some time when the control Pc is zero or near 
zero, the control Ac begins to follow a linear law, as 
suggested in egn. 2.7: 



Ac= Aco 



[’-.V] 



Notice however that the term Aco in this equation is no more 
the actual initial ccmmanded acceleration, but than one the 
missile would have if its initial roll angle was equal to 
the final 

Such behavior defines a boundary in the ccnrrol Ac 
which is clearly seen in 3.29, where the commanded accelera- 
tion is bounded by the curve of the control Ac the missile 
would have if its initial was egual to » tf no 

commanded roll was necessary to reach the rargsr) . 

Althoug the missile commands roll angles larger than 
the ideal initial rcll angle, it is possible to do a 
prediction- with no no computer work- of an aproximation for 
the maximum acceleration rhe missile would experience during 
its fligth, as following (see fig. 3.32): 







1 




(3. 37) 


t 

^ Tl 

[ 


COS ^ 


^icieoi 





where cames from eqn.3.50, and Ac is a straight line 

^ ‘dca/ 

as in fig 3.32. 

Other interesting point is than the missile is 
commanding ro reach roll angles aboun twice of ,dett 
the target is at small accelerations, and when an large 
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accalerations, the missile is makiag a small correction on 
its roll angle, as shown in table I. 

One can see from the figures showing miss disnance, 
that the relative position of missile to target is abour the 
same in all three initial cases. Thus, the new ^,jcal 
each point is the same. As the coniputed an t=0 is 

smaller when the target is at small accelerations, the 
correction has to be larger in order no hit the target, as 
seen in figure 3.30. 

In table I, one can see the final miss distances, 
the final miss distance eg to eg and the final roll angle. 
Such results show that with this digitalized model, good 
results has been obtained. 
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Figure 3.1 Representatioa of the System. 
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3.2 ?ariation of Commands with Initial Roll Angle. 
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Figure 3.3 Variation of Roll Angle - Small Angles- 
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Figure 3.4 Variation of Roll ingle *■ Large Angles. 
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gure 3.5 Ideal Initial Roll Angles. 
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Figure 3.6 Scenarios for Simulation. 
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Figure 3.7 Commanded Acceleration- Case 0. 
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Figure 3.8 Commanded Roll Rate- Case 0. 
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Figure 3,9 Miss Distance in Z Direction- Case 0. 
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Figure 3.10 Miss Distance in I Direction- Case 0. 
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Figure 3.11 
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Figure 3,12 Comnauded Acceleration- Case 1. 
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Figure 3.13 Commanded HoLl Rate- Case 1. 
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Figure 3.14 Miss Distaace- easel. 
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Figure 3.15 Roll Angle- Case 1. 
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Figure 3.16 Commanded Acceleration- Case 2. 
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Figure 3.17 commanded Roll Rate- Case 2. 



67 



mss DISTANCE- METERS 

- 100.0 - 80.0 - 60.0 - 10.0 - 20.0 0.0 20.0 - 10.0 60.0 80.0 100.0 



'1 



2ND CASE 

INITIAL TARGET ACCELERATION- -1. G 
INITIAL TARGET POSITlON-100 M 
SAMPLE PER 100-0. 05 SEC 




I 



Figure 3.18 Hiss Distioce- Case 2. 
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Figurs 3.19 Roll Angls- Case 2. 
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Figure 3.20 Commaaded Acceleration- Case 3. 
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Figure 3,21 Commanded Roll Rata- Case 3. 



71 



r 



C«':,r: 

IN IT I Hi. THKGL": RrXij..[:RnTirjN‘- -T. G 
1N:T]RL TfiRbCT PGSlTiCN-lGC N 
Gfl’IPU: PLPIOD-’O.CS S::c 



I 




Figure 3.22 Hiss Distaace- Case3. 
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Figure 3.23 Roll Angle- Case 3. 



73 



■1 



4TH CASE 

INITIAL TARGET ACCELERATION- -4. G 
INITIAL TARGET POSITION— 600 M 
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Figura 3.24 Commanded Accal eration- Case 4. 
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Figure 3,25 Commanded HoLl Rate- Case 4, 
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Figure 3.26 Hiss Distaoca- Case 4. 
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Figure 3.27 Roll Angle- Case 4. 
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Figure 3.28 Corrected aodel. 
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Figure 3.29 Boundary of the Coaniaaded Acceleration. 
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Ftg* 3.30 a Large target accelerations 




Fig. 3.30 b Small target accelerations 



Figure 3.30 Corrections on the Roll Angle. 
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IV. ASSISTS OF £AINS, SAMPLE RATS AND PITCH ANGLE 
A. ANALISTS OF THE GAINS 

In chapter 3 has been developsd a solution for 
optimal control of a systam as: 

x(k + 1)= A (k) X (k) + E(k) u (k) + E g 

It would ba interssring ~o check if the optimal gains 
reach steady stats, but at the momsit chat the sxter.tion for 
large roll excursions has been introduced, and the system is 
being feed with optimal commands which are varyina each step 
of time, such idea can not be applyed. However one can do 
such check in the model for small roll excursions, which is 
acctually used to compute the optimal commands. 

Doing this, one has the time history of those gains as 
in 4.1 and 4. 2 

One can notice in fig. 4,1 that the gain FG(2,1) associ- 
ated with the effect of gravity has no effect on the 
ccmanded roll rate, and that FG(1,1), as shown in fig. 4.2, 
has a large effect on the commanded acceleration. 
Furthermore, this gain reachs steady state very fast, thus 
it can be assumed that the gain FG will be equal to the 
steady state value during all time. 

From eqn. 3.19, and assuming steady state: 

u(k)= -F(k) x(k)- FG(k) g (4.1) 



and substituting g: 
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u (k)= -F (k) X (k) - C 



(a. 2) 



whare the second term in the right hand side is a constanr, 
and its value is exactly equal to tie value of the commanded 
acceleration necessary at t=0 to correct the gravity fall of 
the missil (or to correct the initial ZSM due to gravity) . 

It might be supposed that one could solve the optimal 
control problem for the system represented by eqn. 4.1, just 
considering one reduced system represented by: 

x(k+1)= A (k) X (k) ♦ B(k) u (k) (4.3) 

with a bias in the ccntrol as: 

u (k)= -F (k) X (k) ^ C (4.4) 

But as shewed in the following analysis, this is not 
po s s i bl 2 . 

The constant term in the right hand side of eqn. 4.4 is 
calculated as follows: 
from eqn.: 

initial ZEM due to gravity= ZEM = " g t^ 



and the gains F(k) are calculated using a Sicatti equation 
as usual. 



C 



Z’t Hg 
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Case 5 was tested using the above considerations, and 
using the same scenario of case 3. 

Figures4,3 andU,4, show the time history of the 
controls. One can see that the commanded acceleration has 
begun at same values as in case 3, but the commanded accel- 
eration reachs a peak considerable higher, them dcreases 
does not following a linear law, with a final at 14.7 meters 
per second square, being this terminal value due to the 
constant term representing the affect of the gravity. 

Refering to fig. 4.4, the commanded roll rate begins at 
a same value as in case 3, but as the control Ac is too 
high, it reachs negative values, going to zero almost at the 
end of the running time. This behaviour of the control leads 
to a large miss distance as seen in figure 4.5, and tablell. 

Figure 4.6 shows the time history of the roll angle, 
which rises to values close to 1.5 radians. As the accelera- 
tion at this point is larger than the correct value, the the 
corrections are excessive and the roll angledecr eases at the 
end of the running time to the value of .12 radians. 

Thus, one can see that the gain due to the gravi-^y's 
acceleration can not be replace! by its steady state 
value. This kind of simplification can thus not be done in 
the system being studied. 



B. EFFECT OF THE SAMPLE RATE 

Throughout all the simulatios a sample period of .05 
seconds has been used. In this section a brief study of the 
effect of the change of this sample rate is given. 

Two best cases have been selected to ilustrate the 
effect of the sample period. 

The first case, case 6, has been run with a sample 
period of .1 seconds and consits of the same scenario as 
case 3. 
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As one can see in figurs4.7 and taolelll, 
noticable change in the commanded acceleration, bur the 
commanded roll rats begins at smaller value than in case 3, 
as is seen in fig. 4. 8 This initial decrease in roll rare, 
leads to a large miss distance in I direction as shown in 
fig. 4.9, and to a small value of roll angle (see fig. 4. 10). 

The second case, case 7, has a smple period of 0.025 
seconds. There is no noticeable change in the time history 
of the control Ac as shown in fig. 4. 11. The commanded roll 
rats begins at a higher value than in case 3 as shows fig. 
4.12, which leads to a final miss distance in Y direction of 
-2.22 meters and in Z direction smaller than case 6 (see 
fig4.13). Figure4.14 shows that the final roll angle is 
increased and the missile cross the target wirh 1.28 radians 
and with a CG-to-CG distance of 2.5 meters. 

In both cases, the miss distance was increased over the 
nominal value obtained with a sampling rate of 0.05 seconds. 
Thus, it would appear that there is an optimal value for the 
sampling rate, which may be conected with the geometry of 
the scenario and with time to go. 

C. EFFECT OF THE INITIAL PITCH ANGLE 

It is important at this step to remenber that througout 
this work has we have been discussing a dimensional model, 
where there is no information on the X coordinate, so it is 
impossible to verify the behavior of the pitch angle. 

In this work, since in all the previous scenarios the 
initial angle 9 was equal to zero, this value has been kept 
as a constant during all time, and considering that whithout 
any information of a third dimension it was not possible to 
correct the time to intercept, this time was also kept 
constant and equal to the nominal value of 5 seconds. 
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Notice “hat this assumption is liicely to be correct if 
one has the horizontal initial distance from target to 
missile compared with the initial vertical distance between 
target and missile large enough in order to have 
small angles. 

The question that rises is, how could this pitch angle 
affect the system if it was net small? 

As seen in fig. 4. 15, the missile velocity in the X 
direction would be: 



COS 9 + g COS0 sinO t 



(4.5) 



which has an affect not only from the pitch angle, repre- 
sented by the cos 0 , but also from the gravity's accelera- 

tion, which will leads to a different time to intercept. 

Considering the same physical scenario as in case 4, but 
changing the initial pitch angle, in order to have the 
missile pointing to the target (sea fig. 4. 16), and keeping 
the missiles velocity of 1090 m/sec in the X direction, one 
has a ccrapletly different geometry of the problem as seen 
from the flight path reference frame. 

With this new situation (see fig. 4. 16) , case 8 has been 
run. Figures 4.17 and 4.13 show that now the missile is 
comanding large positive accelerations, and the roll rate at 
the begining of the flight is too high, going to zero in a 
very small period of time. The miss distance as seen in 
fig.4.19 are increased in the initial part of the fligth and 
as the missile corrects its trejectory it is decreased to 
reach a final CG-to-CG miss distance of 2. 17 m. The roll 
angle, due to the large control Pc is oscilatory in the 
begining and becames constant with a value of .57 radians 
(see fig. 4. 20 and table IV) . 
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Notice that the high values of acceleration 
in some part due to a vertical component of target’s 
velocity, which is seen from the flight path frame as the 
target was manouvering in the 2 direction with constant 
velocity. These large accelerations leads to roll rates too 
large for the physical integrity of the missile. This means 
that although althoug the good results obtained, if compared 
with case 4, they are not practical. 

In order to get rid of the vertical manouever of the 
target, case 9 has been run. In this case, the scenario is 
the same as before with the target also pointing down, with 
the same pitch angle as the missile, and has a X velocity 
equal to the previous case (see fig. 4. 16). 

From fig. 4. 21, one can see that the decrease in the 
control Ac is substantial if compared with case 8, but the 
commanded roll rate is still toe large as shows fig. 4.22 The 
time history of the miss distance is better, resulting in a 
final CG-to-CG distance of abcut 50S of case 8 (see 
fig. 4.23). The roll angle is not oscilatory as seen in 
fig. 4.24 and the missile cross the target with 1.3 radians 
in roll, (see tablel7). 

The results obtained in the two previous cases, suggests 
that the algorithm developed in this work could be readly 



applied 


to air 


-to-surface mis 


siles. In the 


latter, the 


scsnaric 


would 


be favorable t 


o the missile tha 


n in either 


pr9 vioas 


cases , 


since one can 


consider that the 


target could 



be essencially stationary in comparison with the missile 
speed . 

Case 10 has been run with this assumption, and the scen- 
ario as in fig. 4.16. The target'is with zero velocity and 
acceleration, and the missile begins its flight 600 meters 
above it, with the same initial pitch angle as before. 

Fig. 4.25 shows the time history of the commanded accel- 
eration, where one can see that as there is no roll rate to 
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be commanded, the acceleration is following a straighc line 
with very resoanable values. The miss distance is show in 
figure4.27, which shows the final CG-to-CG distance of .31 
meters. 





Based 


in the results 


0 f t 


hese tests. 


one can see that 


ther 


e will 


be some effect 


of 


the pitch 


angle on the miss 


dist 


ance , 


not only due to 


its 


effect on 


the time to inter- 


espt 


, but 


also because at 


the 


moment that there is a pitch 



angle different from zero, even if the rargez is keeping its 
flight level, in the flight path frame a component of the 
target's velocity will show up leading the missile to 
command large accelerations and roll rates. Althoug this 
results harm the performance in an air-to-air missile, in 
the case of air -to-surf ace missiles, when the target has 
been considered with no motion, good results has been 
obtained. 

D. EFFECT OF TIME TO INTERCEPT 

In the simulation of case 1 and 2 in chapter 3, it has 
been observed that when the target was at small accelera- 
tions, the missile did larger corrections on its roll angle, 
with respect to the ideal initial roll angle, than when the 
target was with large accelerations. One can think that must 
be seme kind of compromise between the velocity rate of 
target and missile (which will reflected on nhe time to 
intercept), and the relative position between them, which 
will affect the miss distance. 

In order to do a brief analisys on nhis, case 11 and 12 
has been run. 

In case 11, the scenario of case 2 has been kept, with 
the exception that the missile's velocity was change to 2000 
m/sec., which means that Ti was changed to 2.5 seconds. 
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Figure 4.29 shows that the icceleration is largely 
increased due to the small time required to correct the Z3M, 
and the commanded roll rate is almost twice of case 2 (see 
fig. 4.30). The final miss distance is more than four times 
the value obtained in case 2 , as seen in fig. 4.31 and table 
V. The final roll angle is about half as in case 2 , since 
the projected final position of the target, in the Y direc- 
tion is less chan in case 2 (see fig. 4.32). 

Case 12 was cun with a scenario less favorable to the 
missile, where all the conditions of the previous case was 
kept, except the target position that has been increased to 
200 meters above the missile. 

Now, the commanded acceleration are much larger, with a 
initial Aco of 103 m/sec^, being almost impossible to see 
the difference of the time history of the acceleration from 
one straight line, as shown in fig. 4.33 The commanded roll 
rate is small, about the same as in case 2 (see fig. 4.34). 
The change in the miss distance is noticeable, with a final 
CG-to-CG distance of 4.8 meters as in table 7 , and figure 
4.35. The final roll angle explain the shape of the acceler- 
ation curve, since with the small roll angle as shown in 
fig. 4.36, the system is behaving as for small roll excur- 



sn 0 n s . 








Notice that 


from this analysis 


, one has to realize t 


hat 


there is seme 


kind of envelope 


where the 2-D system 


is 


valid. And in 


order tc determine 


this envelope, one has 


to 


take in account 


not only the geom 


etry of the scenario. 


bat 



also the time to intercept, which is determined not only by 
the relation of velocities of missile and target, but the 
pitch angle too. 
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Figure 4.1 Gains affecting the commanded acceleration. 
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Figure 4.2 Gains Affecting the Commanded Roll Rate. 
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5TH CASE 

INITIAL TAAGET ACCELERATION- ~i. G 
INITIAL TARGET POSITIQN-IOO M 
SAMPLE PERIOD-0.05 SEC 




Figure 4.3 Commanded accaleration-Case 5. 
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ROLL RATE- RDS/SEC 



5TH CR5E 

INITIRL TARGET ACCELERATION- -4. G 
INITIAL TARGET POSITIQN-IQO M 
SAMPLE PERIOD-0.05 SEC 




Figure 4,4 commanded BoLl Sate-Case 5. 
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5TH CASE 

INJTIPL TARGET ACCELERATION- -4. G 
INITIAL TARGET PDSITION-lOO f1 
SAMPLE PERIQD-0.05 SEC 




Figure 4.5 aiss Distance-Case 5. 
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5TH CRSE 

INITIAL TARGET ACCELERATION- -4. G 
INITIAL TARGET POSITIQN-100 M 
SAMPLE PERlDD-0.05 SEC 




Figure 4.6 Roll Angle-Case 5. 
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TABLE II 

Results from Biased Control 
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6TH CASE 

INITIAL TARGET ACCELERATION- -4. G 
INITIAL TARGET POSITION-IOQ M 
SAMPLE PERlOD-0.10 SEC 

a 




Pigura 4.7 Commanded Acceleration-Case 6. 
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Figure 4.8 Commanded Roll Rate-Case 6. 
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6TH CASE 

INITIAL TARGET ACCELERATION- -4. 6 
INITIAL TARGET PCSITION-IOO M 
SAMPLE PERIOD-0.10 SEC 
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Figure 4.9 Miss Distaoca-Case 6. 
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6TH CASE 
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Figure 4.10 Roll Angle-Case 6. 
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7TH CASE 

INITIAL TARGET ACCELERATION- G 
INITIAL TARGET P05ITI0N-1O0 M 
SAMPLE PERIOD-0.025 SEC 




Figurs 4.11 Commanded Acceleration-Case 7. 
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ROLL RATE- RD5/5EC 
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7TH CRSE 

INITIAL TARGET ACCELERATION- -4. G 
INITIAL TARGET P05ITI0N-100 N 
SAMPLE PERIQD-0.025 SEC 
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Figurs 4,12 Commanded aoLl Rate-Case 7. 
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7TH CfiSE 

INITIAL TARGET ACCELERATION- -4. 6 
INITIAL TARGET POSITION-100 M 
SAMPLE PERIOD-0.025 SEC 




Figure 4.13 Miss distaace-Case 7. 
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ROLL ANGLE- RD5 



7TH CASE 

INITIAL TARGET ACCELERATION- -4. G 
INITIAL TARGET POSITION-100 f1 
SAMPLE PERI QD-Q. 025 SEC 




Figure 4.14 Roll Angle-Case 7. 
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TABLE III 

Results for Differents Sample Periods 
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Figurs 4.15 Effsct cf Initial Pitch Angle. 
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Ftgure 4,16a Case 8 




Figure 4.16b Case 9 




Figure 4,16c Case 10 



Figura 4.16 Scenarios With Pitch Angl 
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8TH CRSE 

INITIRL TRRGET RCCELERRTION- G 
INITIRL TARGET POSITION— 600 N 
SAMPLE PERIQD-0.05 SEC 





?igura 4.17 Commandad accaleration-Case 8. 
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ROLL RRTE- ROS/SEC 
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0TH CRSE 

INITIRL TRRGET RCCELERRTION- -4. G 
INITIRL TRRGET POSITION— 600 t1 
SRMPLE PERlOD-0.05 SEC 




Figure 4.18 Commanded Roll Rate-Case 8. 
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Figure U.19 aiss Disrancs-Case 8. 
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Figure 4.20 Roll Augla-Casa 8. 
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INITIPL TPRGET PCCELERRTION- *4. G 
INITIRL TPRGET POSITION— 600 N 
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Commanded Accaleration-Case 9- 
1 1 2 



Figure 4.21 
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9TH CnSE 

INITIPL TRRGET RCCELERRTION- -4. G 
INITIRL TARGET POSITION— 600 H 
5RHPLE PERlOD-0.05 SEC 







Figure 4.22 Commanded Roll Rate-Case 9. 
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Figurs 4.23 Hiss Distaaca-Case 9. 
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Figure 4.24 Roll Angle-Case 9. 
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INITIPL TPRGEIT PCCELERPTION- -4. G 
INITIPL TPRGET POSITION— 600 f1 
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Figure 4.25 Commanded Acceleration-Case 10- 
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IDTH CASE 

INITIAL TARGET ACCELERATION- -1. G 
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Fignra 4,26 Conaanded Roll Rate-Case 10 



lOTH CASE 
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INITIAL TARGET POSITION— 600 M 
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Figure 4.27 Miss Distaace-Case 10 
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Figure u.28 Soil Angle-Case 10. 
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TABLE IV 

Results Using Pitch Angle 
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Figur© 4.29 Coinniaiided Accslsra-ion-Cass 11. 
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Figura 4.30 Commanded Roll Hate-Case 11. 
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Figure 4.31 Hiss Distance-Case 11. 
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Figure 4,32 Roll Angle-Case 11* 
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Figure 4,35 Miss Distaacs-Case 12. 
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ZINAL CONCLaSIOMS ASD COMMENTS 



The scope of the present work was rhe development of an 
optimal digital control to be applied on a bank-to turn 
missile. 

A two dimensional model, as sugested in reference 1, was 
adopted. Afrer the digitalization of the continuous model it 
was necessary to solve a modified Ricatti equation since in 
the stats equation there was a third term representing the 
gravity's effect. The approach that has been adopted is new, 
and although good results were obtained for the scenarios 
considered in this work, is necessary that rhe algorithm be 
further tesred and evaluated in similar problems due to its 
novelty. 

The optimal was solved with an initial restriction to 
small angles. This condition was later relaxed so that large 
roll angles could be analyzed. 

It is dificult to compare the present work with previous 
results since Stallard has indicated a mistake in his orig- 
inal paper, and further works in this area was not found. 

However some comparison with Stallard work is possible. 
The commanded acceleration of the missile are such as to 
correct the ZEM at each point, this agrees with that refer- 
ence. There is a proportional relationship between the 
commanded roll rate and the commanded acceleration , and the 
commanded roll rate is proportional to the defined ^ ,deul s- 
each point, which again agrees with reference 1. 

The algorithm developed in this work requires extensive 
computation at each step, and is clear that some software 
optimization will be needed. The motivation for considering 
the constant steady-state gain due to gravity was to 
decrease this computational burden. This approach however 
resulted in unacceptable miss distance. 
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Another point of investigation that could reflect or. the 
period available to the computer to perform its calculations 
was a change in the sample rats. Trfo different sample rates 
were investigated, both lead to larger errors than the 
nominal period of . 05 seconds. A detailed study on this 
issue is 'left as sugestion for future works, since some 
optimal value of the sampling rate is clearly indicated. 

It is important to keep in mind that the model adopted 
is two dimentinal, while the actual problem is three dimen- 
sional, thus some brief studies were conducted in order to 
check the region of validity of the 2-D. 

In the analysis of the pitch angle, one can see that is 
necessary to have small variations in pitch in order to 
approximate it as a constant. However, at the moment that 
this angle is different from zero, as explained in chapter 
4, it is possible to have in the flight path reference frame 
a target maneuver in the Z direction that will lead to large 
acceleration commands, leading the missile to large miss 
distances, when considering a movable target. When the 
present system was tested tested against fixed targets, the 
results were quite good, this suggests the application of 
the model in air -to-sur face missiles. 

Further investigation were made on the effect of time- 
to-go. As expected, decreased time to go, results in 
increased miss distance. A detailed analysis of more 
complex scenarios is needed in order to properly define the 
effect of time to go. 

Also, it would be interesting to extend the model to 
three dimensions and i.nclude the effects of lags on the 
system in future works. Finally, in appendix A, the 
computer model used in this work is enclosed. Some improve- 
ments in this program can be done, mainly in the data intro- 
duction, and in some optimization of the running time. 
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